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We propose a novel method of measuring the CKM matrix element Vub, via inclusive B decays to 
"wrong sign" charm. This process is mediated by the quark level decay b — > ucs' . When normalized 
to the inclusive semileptonic decay rate, the theoretical expression is very well behaved, with small 
uncertainties from unknown quark masses and no leading renormalon ambiguity. We compute the 
decay rate for this process, including the leading perturbative and nonperturbative corrections. 



The accurate measurement of the Cabibbo-Kobayashi- 
Maskawa parameter V u b remains an important outstand- 
ing problem in B physics. The magnitude of V u b corre- 
sponds to one of the sides of the Unitarity Triangle, the 
angles of which will be constrained or extracted from ex- 
periments to be performed soon by the CDF, CLEO-III, 
BaBar and BELLE collaborations. Measuring \V u b\ with 
precision would provide a vital complementary check on 
the adequacy and consistency of the CKM framework for 
the physics of flavor in the Standard Model jlj . 

The methods which are currently available for probing 
V u b are unfortunately plagued by dependence on phe- 
nomenological models whose uncertainties are difficult to 
quantify reliably. As a result, despite heroic experimen- 
tal efforts present constraints on this parameter are un- 
acceptably weak. The analyses which have been used fall 
into two classes, inclusive decays of the form B — ► X u lv, 
and exclusive transitions such as B — > (%,p)lv. 

The inclusive decay rate has the advantage that it 
can be predicted in the form of a systematic expan- 
sion in powers of 1/mb However, the measure- 
ment of this process is complicated by the overwhelm- 
ing background of inclusive B — > X c lv transitions. At 
present, this background can be suppressed only by con- 
fining oneself to kinematic regions in which only charm- 
less final states can contribute, such as Ei > 2.3 GeV 
or M(X) < 1.9 GeV. Unfortunately, the same operator 
product expansion techniques which allow one to calcu- 
late reliably the total inclusive rate break down when the 
phase space is restricted in this way. Phenomenological 
models must then be used to reconstruct the rate in the 
unobserved kinematic regions, and the model indepen- 
dence of the analysis is lost HQ] . 

Exclusive transitions such as B — * (tt, p)lv are easier to 
study experimentally. On the other hand theoretical pre- 
dictions of exclusive decay channels are polluted by our 
ignorance of the physics of quark hadronization. Various 
approaches to this problem have been proposed (for ex- 
ample, heavy quark symmetry, chiral expansions in the 
soft pion limit, dispersion relations, QCD sum rules, and 
lattice calculations), but in many of these cases signifi- 
cant model dependence remains. It is therefore desirable 



to continue to look for a model independent method for 
the extraction of V u b- 

We propose that V u b may be extracted from the in- 
clusive nonlcptonic transition b — ► ucs', where s' — 
s cos 9c — d sin 9c is the flavor eigenstate which couples 
to c and we take m s = m c i = 0. Here 9c is the Cabibbo 
mixing angle, with sin 6^ ~ 0.2205. This process stands 
out for its theoretical simplicity. Because the transition 
involves four distinct quark flavors, only diagrams pro- 
portional to V u b enter the theoretical expression, with 
no complications from penguins or strong rescattering 
contributions. (For similar reasons, it also has been pro- 
posed recently that V u b may be extracted from exclusive 
processes mediated by this quark transition ||.) More- 
over, as we show below, the ratio of this rate to that 
for b — > civ receives well controlled radiative corrections, 
with no leading infrared renormalon ambiguity. We will 
compute the largest nonperturbative corrections, which 
are spectator-dependent and arise at order We 
warn the reader that our analysis will rely implicitly on 
the use of parton-hadron duality. This assumption is 
common to all extractions of V u b from inclusive B de- 
cays, and while it is not unreasonable to expect it to 
hold in this case, it has not been proven rigorously to do 
so. 

Is such an analysis feasible experimentally? What is 
required is to find "wrong-sign" charmed states in flavor- 
tagged B decays. This is certainly a challenging mea- 
surement. However, some of the necessary techniques 
already exist, as can be seen in Ref. In addition, the 
combination of moving pairs of B mesons and the excel- 
lent vertexing capabilities of the the BaBar and BELLE 
detectors should certainly improve the prospects for the 
measurement of this transition J7J. Of course, the largest 
source of "wrong-sign" charm is from decays mediated 
by the transition b — ► ccs' , which are enhanced roughly 
by a factor of 100 over those which we propose to study. 
Whether such final states can be rejected at this level 
is an experimental question which we are certainly not 
equipped to answer. However, we note that while numer- 
ically the problem would seem to be similar to the quite 
intractable one of digging b — > ulv out of the b — > civ 
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background, it is actually quite different in its details. 
In particular, in the case of nonleptonic decays there is 
no missing neutrino, and perhaps the enormous samples 
of b — > ccs' decays which will be available to BaBar, 
BELLE and CLEO-III will allow one to understand the 
kinematic features of those decays well enough that they 
may be subtracted reliably. In any case, our purpose in 
this Letter is simply to argue that the theoretical situa- 
tion is so attractive that the feasibility of the experiment 
is worth investigating. 

We begin with the transition b — > cX at lowest order. 
The decay is mediated by the local AB = 1 Lagrangian, 
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where P F is the left-handed projector and N c = 3 is the 
number of colors. At tree level, c% = 1 and c 2 = 0; to 
leading logarithmic order, taking a s (mz) = 0.117, the 
Wilson coefficients at the scale m& are given by c\{m b ) ~ 
1.11 and c 2 (m b ) ~ -0.25 @. 

The inclusive decay rate in the channel of interest may 
be written as the imaginary part of the matrix element 
of the forward transition operator T, 

Y(b^cX) = {B\T{b^cX^b)\B) (2) 

= —{B\lm{i [ d 4 xT[tf(x),L(0)})\B). 
m b I J J 

At lowest order, for which the quark level process is & — > 
ucs' , we parameterize the cs' loop by 



Q 2 9^) + B(q 2 )q tl q y ] , (3) 



where is the momentum carried by the cs' pair. Then 
the decay rate may be written as H 



r(6 -> ucs') 



3N c r \v ub \ 2 xi 



l/x 



dz (l/x c — z) 



x [A{z) {2z + l/x c ) + B{z)/x c ] , (4) 



with x c = ml/ml and L = G 2 F ml / 192tt 3 . At tree level, 
with q 2 = zm 2 , the form factors Aq(z) and Bq(z) are 



A (z) 



Xa (z) = |x(l-lA) 2 (l + l/2z), 



B (z) = X bo(z) = x(l-l/zf/z 



(5) 



with x = (c 1 +c 2 /N c ) 2 + 2cl/N c ~ 1.09. The two terms in 
X come respectively from the squares of the color singlet 



and color octet operators in C. Neglecting the octet op- 
erator would induce only a four percent error in X - Inte- 
grating over z and normalizing the result to the semilep- 



tonic decay rate T(b — 
f(x c ) = l — 8x c 



civ) - 
+ 8x1 



r \V cb \ 2 f(x c ), where 
- xt — 12x 2 hix c , 



yields a very simple expression at leading order, 
R = T(b -> cX)/N c T{b -► civ) 

= xKb/V cb \ 2 {l + 0(a s ,l/m 2 )Y 



(6) 
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Note that at this order the ratio R is independent of 
the heavy quark masses. With no corrections included, 
the only distinction between these two processes is in 
the kinematics of the final state products. However, the 
phase space function f(x c ) turns out to be the same in 
the two cases and cancels in the ratio. The only depen- 
dence on x c will be in the higher order corrections to R. 
Hence the effect of uncertainties and ambiguities in the 
definition and extraction of m c , m b and x c will be much 
smaller in R than in the individual rates. 

We now turn to the effects of perturbative QCD. The 
radiative correction to T(b — > civ) comes from the fi- 
nite renormalization of the clI^l current, plus gluon 
bremsstrahlung. It has been computed analytically by 
Nir [0, and the result can be put into the form T(b — > 
civ) = ro\V cb \ 2 f(x c )[l+a bc (x c ) a s /Tr]. The radiative cor- 
rection to T(b — + ucs') is more complicated, and we will 
only consider the corrections to the color singlet transi- 
tion. In this case we can neglect gluons which propagate 
from the bu line to the cs' line, simplifying the calcula- 
tion enormously. As noted above, at tree level the color 
octet transition only contributes four percent to the total 
decay rate, so the error in omitting the radiative correc- 
tion to this tiny contribution should be small. Of course, 
this is a phenomenological approximation, which is not 
formally consistent in terms of an expansion in leading 
logarithms. Here we sacrifice formal consistency in favor 
of identifying and keeping the terms which are numeri- 
cally the largest. A more complete analysis is presented 
in Ref. Q. 

In this framework, then, there remain two classes of ra- 
diative corrections to T(b — > ucs'). The first, a bu , comes 
from finite renormalization of the ml7' x 6l current. This 
correction, which can be extracted from an analysis of 
T(b — > utv) by Czarnecki, Jezabek and Kiihn |l^| , also 
depends on x c through the phase space of the cs' pair. 
The second contribution, a Es ' , comes from the radiative 
corrections to the cs' loop PJl^|. The effect, in the case 
of the color singlet operators, is to add corrections to the 
form factors in H°^ v {q), 
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FIG. 1. Radiative corrections to R as a function of x c . 
We display curves for a £s / (dashed), aiu (dotted), aic (dot- 
ted-dashed), and the total g(x c ) (solid). 

where fi(z) and f2{z) are functions of the charm quark 
mass which have been computed by Voloshin || . We find 
a radiative correction to R of the form 



R = 



x\V ub /V cb \ 2 {l+g(x c )^ + ...}, 



(9) 



with g(x c ) = ab u (x c ) — o- bc (x c ) + o, Ss '(x c )- The function 
g(x c ), along with each individual term, is plotted in Fig. |l| 
for 4.5 GeV < m b < 4.9 GeV, which via the relation m b — 
to c = tub — mo = 3.34 GeV corresponds to the range 
0.06 < x c < 0.12. For the "central value" x c = 0.09, and 
with a s (m b ) = 0.22, one has g(x c )a s /n = 0.21. We see 
that the one loop radiative correction is fairly large. 

In principle, we could also include corrections of order 
a 2 /3o- Such terms, and all those of the form cv"/3g _1 , play 
an important role in the resolution of infrared renormalon 
ambiguities associated with dependence on the pole mass 
TOfc |l4| . In particular, if m b (or equivalently, the Heavy 
Quark Effective Theory (HQET) parameter A) appears 
in an expression at leading order, then the associated 
perturbative series has a renormalon ambiguity at order 
Aqcd/"^, and the term proportional to a 2 ./3o typically is 
large and must be included. In the present case, however, 
the leading dependence m\f(x c ) cancels in the ratio R, 
and therefore the leading infrared renormalon cancels as 
well. As a result, there is no reason to expect two loop 
corrections to be enhanced, even those proportional to 
a 2 /3o- It is another attractive feature of the process b — > 
ucs' that, when normalized to b — > civ, its perturbative 
series has no leading renormalon and is well behaved. 

We now turn to the nonperturbative corrections to R. 
Because R is independent of m b and m c at leading or- 
der, the first corrections appear at order 1/to 2 . They are 
proportional to the HQET parameters Ai and A2, which 
parameterize bound state effects (kinetic energy and hy- 
pcrfine interactions) of the b quark in the initial B me- 
son These corrections were calculated for the pro- 
cess b — > ctv in Ref. . The results of that calculation 
apply to 6 — > ucs 1 with the identifications to c — > and 



to t — > m c ; they apply to b — > civ in the limit m T — > 0. 
However, inspection of Eq. (2.15) of Ref. |ll| reveals that 
the correction is actually symmetric in m c and to t , and 
hence the corrections at order 1 /to 2 cancel exactly in R. 

The leading nonperturbative corrections arise, then, at 
order 1/to|. Two types of contribution appear at this or- 
der. First, there are additional bound state effects, which 
may be written in terms of six new HQET parameters p\ , 
P2 , and T\ , . . . 71 [jlTj . We have no reason to suppose that 
these contributions should be enhanced over their usual 
size, typically at the percent level in total rates. On the 
contrary, we might expect somewhat of a cancellation in 
R as obtained for Ai and A2. We will not include such 
terms here. The second type of effect comes from an- 
nihilation processes such as bu — > cs'. These diagrams 
contribute only to nonleptonic decays, hence only to the 
numerator of R. In addition, they are expected to be en- 
hanced over other terms at order 1 /to| by a phase space 
factor of 167T 2 , corresponding to a two-body instead of 
a three-body final state. Finally, they contribute differ- 
ently to B~ and B° decays, and can be probed by com- 
paring R in the two cases. (Similar processes are respon- 
sible for the lifetime difference t{B~) ^ t(B°) [||-|20).) 

Annihilation diagrams are associated with four quark 
operators of dimension six. Performing the operator 
product expansion to this order, we find a contribution 
to the transition operator for B~ and B° decays, 
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G 2 F m b 
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|K 6 | 2 (l-x c ) 2 {^ c ( Cl + c 2 /iV c ) 2 



x [(1 + 2x c )O u s _p - (1 + x c /2)Ol_ A ] 
+ 2c\ [(1 + 2x c )T%_ P - (1 + x c /2)T^_ A ] 
+ sin 2 ecN^ + cjN.) 2 

x [(1 + 2x c )O d s _ P - (1 + x c /2)O d v _ A ] (10) 

+ 2 sin 2 e c c\ [(1 + 2x c )T d s _ P - (1 + x c /2)T d _ A ] }, 
where we define [M 
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(11) 



v_ A and O q s _ P are color singlet, while 



The operators O 
Ty_ A and Tj_ p are color octet. Also, note that the 
terms mediating B° decays are suppressed by an addi- 
tional factor of sin 2 9 C , because annihilating the initial 
state requires an operator with q = d rather than q = s' . 

The contribution of % p to B decays is obtained by eval- 
uating the forward matrix elements (_B,j|7^p \B q ). The 
matrix elements of the four quark operators depend on 
nonperturbative QCD and have not been computed in 
a controlled way. Although models such as quenched 
lattice QCD and QCD sum rules can yield estimates of 
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their values, reliable information will probably have to 
wait until full unquenched lattice analyses are available. 
It is convenient to parameterize the matrix elements in 
terms of "bag parameters" Bi and q pOl , 



(B q \O q v _ A \B q ) 
{B q \0%_ P \B q ) 
(B q \T^_ A \B q ) 
(B g \T§_ p \B g ) 



-if 2 

— AJE 



B 2 , 



-if 2 



1 f2 _2 , 
lJB a m B a e 2- 



(12) 



In the vacuum insertion ansatz, we have B% = B% = 1 
and ei = £2 = 0, and only the color singlet operators con- 
tribute to the decay. In an expansion in powers of 1/N C , 
one finds Bi = 0(1), while = 0(1/N C ); in this more 
general case, the color octet matrix elements are sup- 
pressed but nonvanishing. For example, one QCD sum 
rules estimate gives t\ ~ —0.15 and e 2 ~ plfl . In terms 
of these parameters, the contributions to R(B~ , B°) = 
x\Vub/V cb \ 2 [1 + g{x c )a s /K + 5 3 {B- ,B )] are given by 



S 3 (B-) = 



n + c 2 /N c ) 2 
2x c )B 2 - (1 + x c /2)B 1 ] 



W) = 



xN c m 2 b f(x c ) 

x[(l 

+ 2c 2 [{l + 2x c )e 2 - (l + x c /2)e 1 }} 

167T 2 / 2 (1 - X c f 

xN c m 2 b f{x c ) 



'■0c{n c {c 2 + Ci /N c ) 2 



x [(1 + 2x c )B 2 - (1 + x c /2)B x ] 
+ 2c\ [(1 + 2x c )t 2 - (1 + x c /2) £l ] }, (13) 

With x c = 0.09 and f B = 200 MeV, wc find 

5 3 {B-) = -0.46Si + 0.52S 2 - 0.018ei + 0.020e 2 , (14) 
S 3 (B°) = -0.00032^! + 0.00037S 2 - 0.017ei + 0.020e 2 . 

We see that the annihilation process is much larger in 
charged than in neutral B decays. Using factorization, 
with B, = 1 and a = 0, we have (S 3 (B~), S 3 (B )) = 
(0.059, 4.2 x 10~ 5 ); taking instead the sum rules estimate 
e 2 = -0.15 gives (S 3 (B~),5 3 (B )) = (0.062,2.6 x 10~ 3 ). 
Since it depends dominantly on the largely unconstrained 
octet matrix element e^, there is huge uncertainty in 
S 3 (B°); on the other hand, for any reasonable values of 
the bag parameters the annihilation contribution to B° 
decay is below the percent level and is safely negligible. 
By contrast, S 3 (B~) is much larger, but since it depends 
mostly on the singlet matrix elements Bi, it has a smaller 
fractional uncertainty. 

To summarize, we have computed the inclusive branch- 
ing fraction for "wrong sign" charm production in B de- 
cays, normalized to inclusive semileptonic decay. The 
leading perturbative correction g{x c )a s /iT is analytically 
calculable, independent of the spectator quark, and ap- 
proximately 20%. The annihilation contribution S 3 is 



smaller but more uncertain; it enters at the 5% level for 
B~ decays and is negligible for B° decays. There is no 
leading renormalon ambiguity in i?, and dependence on 
quark masses appears only in the radiative and nonper- 
turbative corrections themselves. As a result, this pro- 
cess is an excellent place to extract \V u b\ from a theoret- 
ical point of view. This measurement will be challenging 
experimentally, but given the crucial role of V u b in con- 
straining the Unitary Triangle and the present lack of a 
viable alternative, it certainly is worth pursuing. 

After this letter was submitted, we became aware that 
the attractiveness of this process for probing V u b has been 
noted previously [ p2| , although not discussed in any de- 
tail. 
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